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ABSTRACT 

Many of the known extrasolar planets are “hot Jupiters,” giant planets with orbital 
periods of just a few days. We use the observed distribution of hot Jupiters to con¬ 
strain the location of its inner edge in the mass-period diagram. If we assume a slope 
corresponding to the classical Roche limit, then we find that the edge corresponds to a 
separation close to twice the Roche limit, as expected if the planets started on highly 
eccentric orbits that were later circularized. In contrast, any migration scenario would 
predict an inner edge right at the Roche limit, which applies to planets approaching on 
nearly circular orbits. However, the current sample of hot Jupiters is not sufficient to 
provide a precise constraint simultaneously on both the location and slope of the inner 
edge. 

Subject headings: planetary systems: formation — methods: statistical 


1. Introduction 

Early discoveries of hot Jupiters hinted at a pile-up near a 3-day period, but recent transit 
surveys and more sensitive radial velocity observations have discovered planets with even shorter 
periods. The data now suggest that the inner limit for hot Jupiters is not defined by an orbital 
period, but rather by a tidal limit, which depends on both the separation and the planet-star mass 
ratio (Fig. 1). This would arise naturally if the inner edge were related to the Roche limit, the 
critical distance within which a planet would start losing mass (Faber et al. 2005). The Roche 
limit separation, an, is given by Rp = 0.462 o/jwhere Rp is the radius of the planet, and 
fj, = m/M* is the planet-star mass ratio. 

The many formation scenarios proposed for hot Jupiters can be divided into two broad cate¬ 
gories. The first involves slow migration on quasi-circular orbits, perhaps due to interaction with a 
gaseous disk or planetesimal scattering (Murray et al. 1998; Trilling et al. 1998). This would result 
in an inner edge precisely at the Roche limit. The second category invokes tidal circularization 


- 2 - 


of highly eccentric orbits with very small pericenter distances, following planet-planet scattering 
(Rasio Sz Ford 1996; Weidenschilling & Marzari 1996; Ford et al. 2001; Papaloizou & Terquem 2001; 
Marzari & Weidenschilling 2002), secular perturbations from a wide binary companion (Holman et 
al. 1997; Wu & Murray 2003), or tidal capture of free-floating planets (Gaudi 2003). These would 
result in a limiting separation of twice the Roche limit, assuming that circularization can take place 
without significant mass loss from the planet^ (Faber et al. 2005; Gu et al. 2003; Rasio et al. 1996). 


2. Statistical Analysis 


To constrain rigorously the distribution of hot Jupiters, we adopt a Bayesian framework, where 
the model parameters are treated as random variables to be constrained by the actual observations. 
To perform a Bayesian analysis it is necessary to specify both the likelihood (the probability of 
making a certain observation given a particular set of model parameters) and the prior (the a 
priori probability distribution for the model parameters). Let us denote the model parameters 
by 9 and the data by d, so that their joint probability distribution function (PDF) is given by 
p{d, 6) = p{6)p{d\6) = p{d)p{6\d). Note how the joint PDF is expanded in two ways, both expressed 
as the product of a marginalized PDF and a conditional PDF. The prior is given by p{6) and 
the likelihood by p{d\6), while p{d) is the a priori probability for observing the values actually 
measured and p{9\d) is the PDF of primary interest; the a posteriori PDF for the model parameters 
conditioned on the actual observations. The probability of the observations p{d) can be obtained 
by marginalizing over the joint PDF and again expanding the joint density as the product of the 
prior and the likelihood. This leads to Bayes’ theorem, the primary tool for Bayesian inference. 


p{e\d) 


p{d\e)p{e) 

p{d) 


p{d\e)p{e) 
f d9p(dl9)p(9) 
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Often the model parameters contain a quantity of particular interest (the location of the inner 
cutoff for hot Jupiters in our analysis) plus other “nuisance parameters,” which are necessary to 
describe the observations (e.g., the fraction of stars with hot Jupiters in our analysis). Since Bayes’ 
theorem provides a real PDF for the model parameters, we can simply marginalize over the nuisance 
parameters to calculate a marginalized posterior PDF, which will be our basis for making inferences 
about the location of the inner cutoff for hot Jupiters. 


^This is very easy to show: consider a planet on an initially eccentric orbit, with initial eccentricity e and pericenter 
distance Cp. Circularizing this orbit under ideal conditions leads to dissipation of energy but conservation of mass and 
angular momentum. Simply equating the angular momentum of the initial and final orbits gives a final circularized 
radius a = rp(l -|- e) ~ 2rp for e ~ 1. 
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2.0.1. 1-d Model 

We start by presenting a simple 1-d model for the distribution of hot Jupiters. The primary 
question we wish to address is the location of the inner edge of the distribution relative to the Roche 
limit. Therefore, we define x = a/ or, where a is the semimajor axis of the planet and ur is the 
Roche limit. We assume that the actual distribution of x for hot Jupiters is given by a truncated 
power law, 

, xi< X < Xu, (2) 

and zero elsewhere. Here 7 is the power-law index and xi and Xu are the lower and upper limits 
for X. The lower limit, xi, is the model parameter of primary interest, while 7 and Xu are nuisance 
parameters. Therefore, our results are contained in the marginalized posterior PDF for x;. 

For simplicity we restrict our analysis to the subset of known extrasolar planets discovered by 
complete radial-velocity surveys, extremely unlikely to contain any false positives. To obtain such 
a sample, we impose two constraints: P < Pmax, where Pmax is a maximum orbital period, and 
K > Krain, where RTmin is a minimum velocity semi-amplitude. We use iPmin = 30m/s, following 
Gumming (2004). We typically set Pmax = 30 d, even though radial-velocity surveys are likely to be 
complete for even longer periods (provided K > K^\u)- This minimizes the chance of introducing 
biases due to survey incompleteness or possible structure in the observed distribution at larger 
periods. By considering only planets with orbital parameters such that radial- velocity surveys are 
very nearly complete, our analysis does not depend on the velocities of stars for which no planet 
has been detected. Note that our criteria for including a planet may introduce a bias depending 
on the actual mass-period distribution. We will address this with a 2-d model below. Note that, 
in this paper, we exclude any planet discovered via techniques other than radial velocities (e.g., 
transits), even if subsequent radial-velocity observations were obtained to confirm the planet. 

Initially, we make several simplifying assumptions to allow for a simple analytic treatment. We 
assume uniform priors for each of the model parameters, ^( 7 ) C^(7min,7max) andp{xi,Xu) const, 
provided xn < xi < Xu < Xuu and zero otherwise. The lower and upper limits are chosen to be 
sufficiently far removed from regions of high likelihood that these choices do not affect our results. 
We assume that the orbital period (P), velocity semi-amplitude (K), semi-major axis (a), stellar 
mass (M*), and planet mass (msini) are known exactly based on the observations. 

We begin by assuming that sini = 1 (orbital plane seen nearly edge- on) for all systems and 
that all planets have the same radius, Rr. With these assumptions, the posterior probability 
distribution is 

n 

p{xi,Xu,^\xi,...Xn) ^ X])~'^'[\x]~^, (3) 

provided that xu < xi < X(i) < < Xu < Xuu and 7 min < 7 < 7 max- Here n is the number 

of planets included in the analysis, is the smallest value of x among the planets used in the 
analysis, and is the largest value. The normalization can be obtained by integrating over all 


p(x|7,xz,x„)dx = x^ ( — 

X 
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allowed values of xi, and 7 . We show the marginal posterior distributions after integrating over 
the nuisance parameters, Xu and 7 , in Fig. 2 (dotted line), assuming Rp = 1.2 Rj. The distribution 
has a sharp cutoff at and a tail to lower values reflecting the chance that xi < X(i) due to the 
finite sample size. 

Next, we adopt an isotropic distribution of inclinations (cosi ~ [/[—1,1]), but we use the 
measured value for radial-velocity planets where the orbital inclination has been determined via 
transits. We show the marginal posterior distribution for xi in Fig. 2 (solid line). The sharp cutoff 
at is replaced with a more gradual tail, reflecting the chance that sin i < 1 for planets with the 
smallest values of x. 

Now consider the consequences of allowing for a distribution of planetary radii. For transiting 
planets we use a normal distribution based on the published radius value and uncertainty. For non¬ 
transiting planets, we assume a normal distribution with standard deviation app. We show the 
resulting marginalized posterior PDFs in Fig. 2. Allowing for a significant dispersion broadens the 
posterior distribution for xi and results in a slight shift to smaller values. We have also explored the 
effects of varying the model parameter Pmax from 8 d to 60 d. We find that this does not produce 
any discernible difference in the posterior PDF for xp 

Our results are sensitive to the choice of mean radius for the non-transiting planets. In Fig. 3 
(upper) we show the posterior PDFs for various mean radii, assuming app = 0.1 Rj. Since few 
planets have a known inclination, there is a nearly perfect degeneracy between Rp and xp Even 
when we include transiting planets, this degeneracy remains near perfect, i.e., p{xi\Rp,xi, ...,x„) ~ 
p(xi ■ -j^lR'pjXi, ...,Xn). However, it remains extremely unlikely that x; ~ 1 for any reasonable 
planetary radius. 


2.0.2. 2-d Model 


We can improve our analysis by more properly considering the joint PDF in orbital period and 
planet-star mass ratio, which we write as 


f^i -fmin) .Hmax) Omim /^max; c) 


p p 


( 4 ) 


provided prmn < p < //max, P < Pm ax . and a{P,M^,) > xi ■ ap{Rp,p). Here a and /? are new 
power-law indices. Again our results are not sensitive to the nuisance parameters, //min, //max, and 
Pmax- For definiteness, we fix their values at Pmax = 30 d, //min = 3.3 x 10 “^, and //max = 0 . 01 . We 
take priors uniform in tan“^(a) and tan“^(/3), as the density P[—7r/2,7r/2] corresponds to uniform 
prior density for the slope of the power-law distribution on a log-log plot. We take a prior uniform 
in logc, as is standard for scale parameters. Our calculation of the likelihood is similar to that 
of Tabachnik & Tremaine (2002), except we replaced their inner boundary of P > Pmin with our 
boundary a(P, M*) > xiap{Rp, p). The necessary integrals can be performed analytically provided 
we approximate 1 -|- // ~ 1. For convenience we used Maple to obtain an analytic expression for the 
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integral of the likelihood over P and The remaining integrals over a, f5, and c were performed 
numerically over a grid with 10^° points. 

By considering the joint mass-period PDF, we are able to account for the bias previously 
introduced by imposing K > iPmax- Since we are considering only planets with orbital parameters 
such that radial-velocity surveys are very nearly complete, our results depend only on the number 
of surveyed stars for which no planet has been discovered, and not the observed velocities of these 
stars. For the total number of stars in radial-velocity surveys that are complete for K > i^min and 
P ^ T’max we estimate = 2000. We show the resulting marginalized posterior PDF for xi in 
Fig. 3 (lower). 


2.0.3. The Shape of the Inner Cutoff 

The above results clearly demonstrate that the present observations strongly favor an inner 
cutoff at the ideal circularization distance rather than at the Roche limit, assuming that the inner 
edge follows the slope of the Roche limit. We have also performed calculations treating this slope 
as an unknown model parameter. Unfortunately, the present observations are not sufficient to 
constrain this parameter empirically, and the resulting marginalized posterior PDF for xi still 
allows, but no longer exclusively favors, xi ~ 2. 

We could gain additional leverage by including planets in short-period orbits down to lower 
Rmin- Unfortunately, the incompleteness of the radial-velocity surveys is likely to be significant for 
Rmin < 20 m/s. Since the number, quality, and spacing of observations varies widely among the 
stars in these surveys, it would be necessary to calculate the probability for detecting planets as a 
function of orbital period and velocity semi-amplitude for each star. Additionally, the planet mass- 
radius relation, which is extremely flat for planets near IMj, becomes important for much lower- 
mass planets. Therefore, we have not attempted to extend our analysis to planets for which radial- 
velocity surveys are not yet complete. We expect that the recent improvements in measurement 
precision will eventually extend their completeness to smaller iFmin- 

We have begun a preliminary investigation of the constraints obtained by adding information 
from the OGLE transit survey (Udalski 2002). Due to both signal-to-noise and aliasing issues, 
the OGLE survey does not provide a complete sample of short-period planets for a significant 
fraction of parameter space. Therefore, it is necessary to calculate the probability of detecting 
planets with various orbital periods and radii. We estimate these probabilities by taking the actual 
observation times and uncertainties for the 62 transit candidates from the 2002 OGLE observations 
of the Carina field and applying the detection criteria from Pont et al. (2005). While the transiting 
planets are consistent with our above findings, they do not provide sufficient additional information 
to constrain the slope of the inner edge. We look forward to both radial-velocity and transit surveys 
detecting additional lower-mass objects in short-period orbits, so that we may eventually constrain 
the slope of the cutoff empirically. 
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3. Discussion 

The current distribution of hot Jupiters shows a cutoff that is a function of orbital period 
and planet mass. Under the assumption that the slope of this cutoff follows the Roche limit, our 
Bayesian analysis solidly rejects the hypothesis that the cutoff occurs inside or at the present Roche 
limit. This is in constrast to what would be expected if these planets had slowly migrated inwards 
on quasi-circular orbits and with radii close to the presently measured values around 1.2 Rj. If 
confirmed by future analyses of a more extensive data set, this result would be highly significant, 
as it would eliminate a broad class of popular migration scenarios for the formation of hot Jupiters. 

Instead, our analysis shows that this cutoff occurs at a distance nearly twice that of the 
Roche limit, as expected if the planets had been circularized from a highly eccentric orbit. These 
findings suggest that hot Jupiters may have formed via planet-planet scattering (Rasio & Ford 
1996), tidal capture of free floating planets (Gaudi 2003), or secular perturbations from a highly 
inclined binary companion (Holman et al. 1997). Regardless of the exact mechanism, our model 
would require that the hot Jupiters all started on highly eccentric orbits and survived the strong 
tidal dissipation needed to circularize their orbits. A few caveats are worth mentioning. Our 
study addresses the statistical properties of the population of hot-Jupiters and does not attempt to 
advance the state of knowledge of any specific planet. In particular, we adopt average properties of 
an assumed distribution that is analogous to—and derived from—the presently known distribution 
of hot Jupiters, but we do not consider or solve for the specific properties of any individual extrasolar 
planet. Moreover, strongly non-random or non-gaussian effects would be poorly modeled with the 
technique developed here. 

An alternative explanation is that the planets migrated inwards at an early time and arrived at 
their Roche limit on a quasi-circular orbit when their radii were still > 2Rj (Burrows et al. 2000). 
The dissipation process causing the migration must then have stopped immediately afterwards to 
avoid further decay of the orbit as the planets continued to cool and contract. We find this scenario 
unattractive, especially since there is no natural explanation for the factor of 2 in this case. 

Yet another alternative is that short-period giant planets are destroyed by some process before 
they reach the Roche limit. HST observations of HD 209458 indicate absorption by matter presently 
beyond the Roche lobe of the planet and have been interpreted as evidence for a wind leaving the 
planet and powered by stellar irradiation (Vidal-Madjar et al. 2003, 2004). Further theoretical 
work will help determine under what conditions these processes can cause significant mass loss 
(e.g., Hubbard et al. 2005) and whether complete destruction could occur rather suddenly when 
the orbital radius decreases below 2aR. 

Future planet discoveries will either tighten the constraints on the model parameters or provide 
evidence for the existence of planets definitely closer than twice the Roche limit. Additionally, future 
discoveries of transiting hot-Jupiters around young stars could help discriminate between the above 
alternatives. Moreover, new detections of lower-mass planets with very short periods could help 
better constrain the shape of the inner cutoff as a function of mass. In the future, an improved 



statistical analysis could also include such low- mass planets, where surveys are not yet complete. 
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Fig. 1.— Minimum mass ratio vs orbital period for the current observed sample. Planets discovered 
by radial velocity surveys are shown as triangles with arrows indicating 1 — a uncertainties in mass 
due to unknown inclination. The magenta squares were discovered by radial velocity observations 
and have inclinations and radii measured via transits. The blue squares show planets discovered by 
transit searches. The green lines show the minimum mass corresponding to various velocity semi¬ 
amplitudes and roughly indicate where radial velocity surveys are nearly complete (> 30m/s), have 
significant sensitivity (> 10m/s), and are only beginning to detect planets (> 3m/s). The two 
red lines show the location of the Roche limit (aji) and the ideal circularization radius (acUc) for 
a planet with a radius Rp = 1.2 Rj. The red lines do not apply to the lowest mass planets that 
likely have a radius significantly less than 1.2 Rj given their different internal structure. 
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Fig. 2.— Marginalized posterior probability distribution for xi, the lower cutoff for the ratio of 
a planet’s semimajor axis to the Roche limit. Here we show multiple posterior distributions for 
various simplified models. The dotted curve assumes Rp = 1-2 Rj and sini = 1 for all planets, 
while the remaining lines assume the observed inclinations and radii for transiting planets and an 
isotropic distribution of inclinations and mean radius Rp = 1.2 Rj for the remaining planets. The 
widths of the distributions in radii are O.ORj (dotted and solid curves), 0.05i?j (long dashes), 0.1i?j 
dots-dashes), 0.2Rj (short dashes). 










